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Abstract 

CP violation leads to a difference between the parameters g~^ and g~ describing 
the energy distributions of the charged pions produced in the —>■ 
and K~ —>• decays. We study the difference {g^ — g~) as a function of 

the relative contributions from the QCD-penguin and the electroweak-penguin 
diagrams. We hnd that the combination of these contributions in {g^ — g~) is 
very similar to the corresponding one dehning the parameter s' in the ^ 2ti 
decays. This observation allows a determination of the value of {g^ — g~), using 
data on e'. 
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1 Introduction 


Since 1989 it has been known that the magnitude of direct CP violation in the 
Kl —> 2tt decays crucially depends on the relative strengths of the imaginary parts 
of the QCD-penguin (QCDP) and the electroweak-penguin (EWP) contributions 
to the amplitude pp. The reason for this sensitivity is that the contributions to 
s' from the two diagrams have opposite signs and partially cancel one another. 

As the dynamical structures of the amplitudes for —>■ (Stt)^ differ from 

those for > 27r, there is no immediate relation between the strengths of 

direct CP violation in the Kl ^ 27t and the —>■ (Svr)^ decays. In particular, 

it was observed in refs. Hlal that in contrast to the situation in Kl 2n, the 
CP violating effect in —>■ produced by the QCDP contribution is 

enhanced by the EWP contribution. 

However, in the present note, we shall demonstrate that the —>■ 

decays are similar to the Kl —>■ 27i decay in that the EWP contribution cancels 
part of the QCDP contribution. Due to this circumstance, we suggest that a 
simultaneous study of the decays —>■ vr^vr'^Tr^ and —>■ Tr^Tr^vr^ could throw 

new light on the relative strengths of the QCDP and the EWP mechanisms in 
direct CP violation. 

We shall estimate, in the framework of the Standard Model, the CP violating 
contributions to the slope parameters and g~ characterizing the charged pion 
energy distributions in the decays (formerly t' decay). The slope 

parameters are dehned by the expansion 

\M{K^{k) 7r°(pi)7r°(p2)7r=^(p3))r ocl + g^Y + ..., (1.1) 

where 

F = (s3 - So)/m^, Si = {k-pif, SQ = m\/?, + ml. (1.2) 

Our tenet is corroborated by a calculation of the amplitudes to leading non¬ 
vanishing order in a momentum expansion. As was previously found for the 
—>■ TT^TT^TT^ decays w, higher-order corrections do not considerably change 
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the conclusion concerning the relative magnitudes of the QCDP and EWP con¬ 
tributions to the difference {g~^ — g~)T- The role of higher-order corrections in 
the r' decays will be considered elsewhere. 

2 The amplitude 

Our starting point is the AS = 1 effective non-leptonic Lagrangian proposed in 

ref. jl] 

L{AS = 1) = V2Gf sin 9c cos ^ ^ CjOj , (2.3) 

i 

where Oi _6 are effective four-quark operators represented by the operator prod¬ 
ucts 


Oi =SL^f,dL ■ ulIi^ul - sl^^ul ■ UL'^y^ldL ; ({8/}, AJ = 1/2), (2.4) 

O2 =SLlf,dL ■ ■ ULli^dL + ■ d.L'l^idL 

+ 2sL'i^,dL ■ ; ({8^}, AJ =1/2), (2.5) 

O 3 =SL'lfidL ■ UlI^Ul + SLli^UL ' Uilf^di + 2sLl^,dL ■ dLl^idi 

- ZsLi^di ■ sl^^sl ; ({27}, AJ = 1/2), (2.6) 

Oa =SLliidL ■ ULlfiUL + Sl'^^Ul ' UL'^^di 

- SL 7 ^dL ■ dLi^dL ; ({27}, AJ = 3/2), (2.7) 

O 5 =SL'yfj,X‘^dL{ ^ dRl 11 ^°' Qr) ({8}, AJ = 1/2), (2.8) 

g=ii,d,s 

Oq =SLlf,dL{ ^ dRlf^dR) ; ({8}, AJ = 1/2) . (2.9) 

q=Uyd,s 


Among these operators, only O 4 generates AJ = 3/2 transitions. The operators 
originate from the QCDP diagrams. To calculate CP-violating effects, also 
the operators generated by the EWP diagrams must be added, 

O 7 =^s7m(1 + 75)c?( X] egg7^(l-75)g); (AJ = 1/2,3/2), (2.10) 

q=u,d,s 

08 =-12 e,{sLqR){qRdLy, e,= (|-^,-^); (AJ = 1/2, 3/2). (2.11) 

q=u,d,s 
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The corresponding Wilson coefficients are small, being proportional to Oem- 
The coefficients c^s contain the imaginary parts necessary for CP violation. 
Bosonization of the operators Oi is achieved through the substitutions [S] 

<?i(l + 75)qk = - ’ (2-12) 

T" F 

+ 75)gfc =i[{d^U)U^ -U{d^U^) - - {d^U)m)]kj . (2.13) 

Here, m is the diagonal quark-mass matrix. 


m = Diag{m„, m^, mJ , 
and the remaining parameters are dehned as 

r = 2ml/{mu + rrid), A 1 GeV, = 93 MeV. 


The 3x3 G-matrix is written as an expansion 

TT A i\/27f f iTc f \ 

7i( ^ A'F) 


in the pseudoscalar nonet-meson-held matrix tt 


71 7! 7! 


TT = 


TT 


K+ \ 


TT 


K- 


ZA _L ZA _ ^ /yo 
73 x /6 72 

-0 TTq 2 'Kg 

73 76 / 


+ 


(2.14) 


(2.15) 


The PCAC condition demands 03 = 0 [HI and we adopt this condition as well, 
bearing in mind, that on mass shell, the values of the mesonic amplitudes are 
independent of the parameter 03 . 

In the calculation of the iC —> dvr amplitudes we make use of the Fierz iden¬ 
tities for the colour matrices 


\o\7 _16X0X7 l\o\7 
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as well as the Fierz identities for the Dirac matrices 


S 7 m (1 + l5)d ■ - 75)g = -2s(l - 75)? ■ q{l + '^^)d . 


Thus, in leading order non-vanishing approximation our result for the matrix 
element can be expressed as 


M ^ 7r°(pi)7r°(p2)vr=^(P3)) = 


= K 


1 

1 ± %aKM H- — 


m 


K 


9C4 


1 _ ^ ) Y{l±ibKM) + ••• 

/Co 


(2.16) 


The kinematic variable Y is defined in eq. (HI. The overall strength is regulated 
by the parameter 

(2.17) 


GFm\ . 

K = -^ Co smb'c; costae , 


6^2 

and the remaining parameters are functions of the following combinations 

32 

Co = Cl - C2 - C3 - C4 7 - /3Rec5 , 

y 


(2.18) 


C 5 — C 5 + 777^6, £7 — C 7 -f- SCg 

16 


(2.19) 


P = 


2 mt 


( 2 . 20 ) 


A2(m„ -7 may 

The terms axM and bxM in the amplitude (jnsD are the imaginary parts of 
the amplitude generated by the Kobayashi-Maskawa phase. Explicitely, we find 


axM 

bxM 


(d 

(d 


/32^ ~ 
Imcs + 

/32^ ~ 
ImC5 + 


6 A^Imc 7 \ 

2 ) 

3A^Imc7 \ , 
772 77T ) / 



( 2 . 21 ) 

( 2 . 22 ) 


with coefficients as above. 

Our approach can also be used to caculate the iF —^ 27r amplitudes. For their 
real parts we get 


M{K\ TT+TT-) 
M(A:+ ^ 7r+7r°) 


GfFtt 

GpF-K 


sin 0(7 cos 9c {my 
sin 9c cos 9c {my 


ml)co , 


ml){lcy. 


(2.23) 

(2.24) 
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A comparison between the real parts of the amplitudes of eqs (irrai and Trm 
shows that their ratio is nothing more than a reflection of the well-known relation 


M{K+{k) vr°(pi)7r°(p2)vr+(p3)) = ^ ) 


X 


1 + 


3mz 


m 


K 


1 - 


3M{K+ 7r+7r°) 

M{Kf TT+TT-) 


Y 


(2.25) 


obtained earlier using soft-pion techniques and current algebra. 

From the data on the K ^ 2n decay rates nm, it follows that 

32 


Cl — C2 — C3 
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■/?Rec 5 = -10.13 , 
C 4 = 0.328 . 


(2.26) 

(2.27) 


Furthermore, the combination /3Rec5 can be determined separately, provided we 
are willing to accept the estimate of Shifman et al. mim, 

Cl — C 2 — C 3 = —2.89 , (2.28) 

which leads to the value 

32 

—/3Rec5 =-7.24 . (2.29) 
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To estimate CP-odd effects in (3vr)^ decays, described by the coeffi¬ 

cients axM and by bKM, we need a certain combination of the parameters Imcs 
and Imc 7 . The theoretical preditions for these parameters are very uncertain and 
different authors (see [3]) give different results. Fortunately, the combination en¬ 
tering the —>■ TT^TT^TT^ amplitude turns out to be similar to the combination 

determining the parameter e' in 77^ —> 27r decay |2]. This circumstance allows us 
to obtain a reliable estimate of — g~)r'- 


3 Estimate of the CP-odd difference [g^ — g )^i 

Although the amplitude (ITT^ incorporates the imaginary terms necessary for 

CP violation, this is not sufficient for producing observable CP-violating effects. 

^In 1^ there is a misprint: a factor y is missing after the factor (1 -|- j^) in the expression 
for the K'^ 7r°7r°7r+ amplitude, eq.(6.9). 
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In fact, the observable effects arise from the interference between these CP-odd 
terms and the CP-even imaginary terms created by the strong-interaction hnal- 
state rescattering between the pions. The strong-interaction effects are intro¬ 
duced into the K ^ Stt amplitudes of eq. (iTlel) by adding two terms, Qr' and br', 
so that 


M{K^{k) 7r°(pi)7r°(p2)vr=^(P3)) 




(3.30) 

F (1 + ± tibKM - aKu)) + ...]. 

V 2co/ 

This assumption is valid as long as the rescattering contribution can be treated 
in the linear approximation. 

The slope parameters and g~ were dehned in eq. (HH). From this dehnition 
and eq. (HOnil we get for the relative difference in slope parameter for the iC —>■ dvr 
decays 


g~^ — g~ \ _ dr'jbKM — 0 >Km)t' 

g^ + g~ ) r' 1 + ar'br' 

The strong-interaction-rescattering parameters, a,-' and br'-, are determined by 
calculating the imaginary parts of the loop diagrams of Fig. 1. Putting the 
intermediate pions on shell (see Appendix) yields, in leading approximation. 




ttr' = 0.12, br' = 0.49. 


(3.32) 


The CP-odd numerator of eq. (IH3TD can be calculated from the expressions in eqs 
and Trm . and is found being equal to 


{bKM — aKM)T' — 


16 c 4 /d Imcs 6 /dA^ Imcy 
Co(Co - |C4) 


m|-co 


1 - 


com^ 


2{m 


m: 


Fco - |C4) 


= 0.042/9 Imcs (1-|-27.8 Imcr/Imcs) . (3.33) 

The combination of Wilson coefficients in this formula, 

f3 Imcs (1 -|- 27.8 Imcy/Imcs) , (3.34) 
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is very similar to another combination 


1-4 


, 24.36 Imdj, (1.63 ± 0.25) ■ 10“ 

/3Imc,(l + —.—) =-- PRec, (3.35) 

defining the direct CP-violating parameter e' in 2tt decay [21 El- The 

parameter takes into acconnt isospin-breaking contribntions generated by the 
two-step transition —>■ 7T^ri{T]') 

At 12 = 0.124 expressions and (HOHll coincide, giving 

Ag^, = ( 1.8 ±0.28) ■ 10 “® , 


(3.36) 


and at 12=0.25 


Agr ' = 2.1 ■ 10“^(1 


4.7 Imcy/Imcs 


1 ± 32.48 Imc 7 /Imc 5 ' 

Both valnes of 12 are in line with estimates fignrating in the literatnre (see 
and references therein). 


(3.37) 


4 The CP-odd difference — g 


As we shall now show, onr resnlt for the slope-parameter difference in t' decay, 
as embodied in eq. (ESI, enables ns to draw qnite precise conclnsions concerning 
the magnitnde of the slope-parameter difference in another decay, namely the 
—>• TT^TT^TT^ decay, or r decay. From eqs (33) and (34) in ref. [Sj, we derive 
the following relation 


{^KM — 0'Km)t — —2 


I 6 C 4/3 Imcs ^ 3/3A^ Imcy 


mj^co 


1± 


12c4m|- 
A^(co ± 9 C 4 ) 


(4.38) 


.Co(co ± 9C4) 

The slope-parameter difference Agj. is again given by expression (ESH), provided 
index t' is everywhere replaced by r. The valne of the rescattering parameter a 
does not change, bnt that of h does. 


ar = 0.12, br = 0.714. 


(4.39) 


Combining eqs (ESH), dSSSl) and dOHD we can form the ratio of the slope 
parameters differences. 


— Agr 2^0 — 904/2 1 — 14.34 Imc7/Imc5 1 + ar'hr' 
Agr' Co ± 9 c 4 1 ± 27.8 Imc 7 /Imc 5 1 ± Urbr 


(4.40) 
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Now, only negative values of the ratio Imcy/Imcs appear in the literature j2]. If 
furthermore, we assume that the numerical value of this ratio is so small that the 
sign of the right hand side of eq. (loni) is positive, we may conclude that 


> 0.56 ■ 10"^ (4.41) 

This result differs from other estimates, as exemplihed by refs. HB] and m 

-Agr = l.SAgr^ |IS|, -Agr = 2.2Agr' (Hj. (4.42) 

Moreover, as follows from our discussion in Sect. 3, we strongly believe Ag^-i 
to be of order 10“®. In contrast, Ag^ can reach values of order 10“^ , providing 
the EWP contribution cancels out a considerable part of the QCDP contribution 
(see eq. onni)). For example, if the EWP cancels half of the QCDP contribution, 
then 

-Ag^ = 7.8Ag^, > 1.4 ■ lO”® (4.43) 

and if the EWP cancels three-quarters of the QCDP contribution, then 

-Ag^ = 17.2A^^/ > 3.1 ■ 10"^ (4.44) 


Ag-r (in units 10 ^) 

Agr^ (in units 10 ^) 

Refs. 

-700 ± 500 

-15 ± 275 

[151 

\^9t\lo < 0.7 

- 

[16] 

-0.16 

- 

[17] 

\Ag,\ = 38.2 

|Af7.,| =31.5 

[18] 

-0.23 ±0.06 

0.13 ±0.04 

[13] 

(—4.9 ± 0.9) sin 5 

- 

[2] 

-2.4 ± 1.2 

1.1 ±0.7 

[14] 

-(3.0±0.5)x; 0.5 < a; <5.0 

- 

[3] 

(-Af7Qio>(0.56±0.09)/(a;) 
At X = 1, (—A^f,-) = 2.9 ± 0.6 

0.18 ±0.03 

present 

work 


Table 1: Values for the slope-parameter ratios Ag^. and Ag^-i in r and t' decays, 
in units of 10“^. 



















These examples show that a simultaneous measurement of Ag^-' and Ag^- 
can clear up the question about the true relative strength of EWP and QCDP 
mechanisms in direct CP violation. The estimates of the values Agr and Ag-r' as 
obtained in other investigations are summarised in Table 1. 

5 Concluding remarks 

We have calculated the CP-odd difference of slope parameters, Ag^f of eq. (ESH), 
in the t' decays —>■ in leading non-vanishing approximation in a 

momentum expansion of the decay amplitude. 

We observe that the difference of slope parameters Ag^-i in 
decay and the parameter e' in —>■ 27r decay both depend practically on one and 

the same combination of the Wilson coefficients Imcs and Imcr. This observation 
permits a reliable estimate of Ag^i using the known magnitude of s'. 

A comparison with the value of the corresponding parameter Ag^ in the r 
decays —>■ shows that Ag^- is expected to be at least 3 times larger 

than Ag'^. In fact, it may be even one order of magnitude larger, provided there 
is a sizeable cancellation between the electroweak-penguin and the QCD-penguin 
contributions to the parameter e'. Such a cancellation is not excluded pen]. 

We have not considered the possibility of a sequential decay —>• —>■ 

TT^TT^TT^ through an intermediate g ^ transition, a correction which is of order 
We shall study this possibility elsewhere. In the case of 
decay, higher-order corrections increase Ag-r by 20%, but change very little the 
relation between electroweak-penguin and QCD-penguin contributions |Hj. We 
expect a similar increase of Ag^-i in the —>■ decay, since ar> ~ Ur, and 

according to ref. [3] corrections increase the value of a,- by 30%. 

Acknowledgments. We would like to thank the Swedish Research 
Council for financial support. One of us (E.Sh) would also like to acknowledge a 
partial financial support from the Grant RFBR-02-02 16957. 
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6 Appendix 


Here, we shall calculate the CP-even imaginary part coming from the pion- 
rescattering diagrams displayed in £g. 1. The imaginary part of a diagram is 
obtained by cutting the internal lines as shown. 


7r°{pi] 




■^°{PuP2) 




■n+{p3. 


a) 


b) 


c)+d) 


Figure 1: Rescattering diagrams for the imaginary part. Diagrams are cut along 
the dashed line. Diagrams c) and d) are related through 7 r°(pi) 7 r°(p 2 )- 

We start from the amplitudes 

M ^ 7r°(pi)7r°(p2)7r+(p3)) = A + B{so-s^) (6.45) 

M (R'+(/c) ^ 7r+(pi)7r+(p2)vr"(p3)) = M + B'{so-s^) (6.46) 

with kinematic variables as defined in eq. (HI. The t' decay amplitude, eq. (16.451). 
is given in eq. dnni), of which we only need the leading real term. In the r 
decay amplitude, eq. ()6.46|) . the parameter A' is twice as large as A. For the tttt 
scattering amplitudes we insert the leading-order approximations. 


M (7r+(gi)7r“(g2) 

7r°(g3)7r°(g4)) = 

^ TV 

(6.47) 

M (vr°(gi)7r°(g2) 

^ 7r°(g3)7r°(g4)) = 

TT 

(6.48) 

M (7r°(gi)7r+(g2) - 

7r°(g3)7r+(g4)) = 

Ut-ml) 

■^TT 

(6.49) 

(?i + ? 2 )^ t = {qi 

- qsf and u = {qi - 

^ 4 )^ as usual. 



First, we calculate diagram a) with a tt+tt pair in the loop. The result is an 
imaginary contribution to the t' decay 


5Ma = 




(s 3 — ^^)\l 1- — \A' + ^B\s 3 — So)] 

S 3 


(6.50) 


10 








However, we are not interested in the exact value of The slope parameters, 
eq. (HH), are dehned through an expansion in H = (ss — so)/m‘l. Moreover, we 
normalise the —>■ tt+tt+tt” decay parameters as 


A' = 


2 2 

^rriK 


B' = 1 + 9c4/co = 0.718 , 


(6.51) 

(6.52) 


so that a short algebraic calculation gives as result 


= i 


. /So - 4/^2 


727rF2 


So 


1 + 


3(s 3 — So) 


m 


K 


1 + 


3so(so - 4m2) 4 


(6.53) 


Diagram b) with two neutral pions in the loop give a contribution to the 
imaginary part 


The parameters for the decay are 


A = \m\ 


B = -(1-9c4/2co) = -1.14 . 


(6.54) 


(6.55) 

(6.56) 


The expansion of this contribution yields the result 

S 3 — So 


5Mb = i 


mj^m^ /so-4/x2 


967rF2 


So 


1 + 


mi 




3B 


(6.57) 


OR i^3so(so - 4m2) 

There are two contributions, diagrams c) and d), with in the loop, since 
the hnal state is symmetric in the two neutral pions, 7r‘^(pi) and 7r°(p2)- We shall 
not give the exact expressions, corresponding to eqs ()6.50|) and ()6.54D . since they 
are somewhat complicated. The expansion of the sum of the two amplitudes 
results in an imaginary contribution 


5Mc + 5Md =i 


m% 


'so - 4/^2 


1447rFf \i So 
3(s 3 - So) 


-( 1 - 


3mt 


m 


K 


+ 


2m\ 


^ ^ 2 m^(so - 2ml) ^ 3m^ ^ 


so(so - 4m2) m 


K 


(6.58) 
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